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Abstract 

The thermodynamic hmit of the Lipkin model is investigated. While the limit turns out to be 
rather elusive, the analysis gives strong indications that the limit yields two analytically dissociated 
operators, one for the normal and one for the deformed phase. While the Lipkin Hamiltonian is 
hermitian and has a second order phase transition in finite dimensions (finite particle number), 
both properties seem to be destroyed in the thermodynamic limit. 
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INTRODUCTION 

The Lipkin model Q], even though it appeared more than forty years ago in the hterature, 
has, up to present, given rise to a vast amount of research papers. This continuing interest 
is due to its basic properties in that it demonstrates the mechanism of a phase transition 
including symmetry breaking for a schematic many body system. Originally introduced as 
a model in nuclear physics [2] it has found applications in a broad range of other topics: 
statistical mechanics of quantum spin systems |3|], Bose-Einstein condensates ^ as well as 
quantum entanglement |5i, to name but a few. Recent interest is focused upon its large 
behaviour with being the particle number. Progress has been achieved using a variety of 

n Q 

methods such as continuous unitary transformations ^ , a semiclassical approach |7| or the 
Bethe ansatz ^. 

We briefly list here the major results relevant to the present paper. In its original form 
the Lipkin model considers interacting Fermions occupying two fi-fold degenerate levels. 
One major appeal is due to its easy solubility. This is achieved by rewriting the kinetic 
energy term of the Fermion Hamiltonian, viz. 



k.a 



and the two body interaction term 

Et t 

k,k',a 

in terms of the SU(2) operators Jz and J^, respectively. In dimensionless form the Hamil- 
tonian thus reduces to 

m) = j^ + ^{Jl + j-) (1) 

where the 2j + 1 = + 1-dimensional representations of the SU(2) operators are used 
with N < Q being the particle number. Here the interaction is scaled by to ensure that 
H is extensive, the operators and effectively scale as A^^. The Hamiltonian allows 
reduction into two spaces: m integer and m half-integer, with m the eigenvalues of J^; 
it corresponds to A^ even and odd, respectively and is denoted as parity. For A > the 
A^-even and A^-odd levels are obviously separated and remain so for all A < 1 while the 
levels become degenerate for A > 1. The phase at A < 1 is called the normal phase while 
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the symmetry (parity) breaking phase at A > 1 is called the deformed phase. For A > 1 
this phase transition is confined to the lower part of the spectrum in that higher up in the 
spectrum the normal phase still prevails. In fact, for a specific A > 1 there is an energy 
i?fc^(A) such that the ordered spectrum is for k < associated with the deformed phase 
while for k > kc the energies relate to the normal phase |9|. In other words, the transition 
moves up in the spectrum from the ground state Ei to En/2 when A increases from unity 
to infinity (we confine ourselves to the lower half of the spectrum as the symmetry of the 
model yields the same pattern in the upper half but mirror refiected). For finite A^, Ej^^ is 
clearly characterised by the minimum of the level distance; this minimum gap vanishes for 
iV — > cxD. In fact, at A = 1 the large N scaling behaviour p, is found 

^E,~U) (2) 



whereas for A > 1 the gap around E^^ behaves as 



We mention that a mean field approach yields the result 



Ek{X) = kVT^ for A<1 (4) 



EkiX) = k v/2(A2 - 1) for A > 1 (5) 

being obviously unable to reproduce the richer structure described above; however, the last 
two equations do describe the lower part of the spectrum remarkably well. 

In the present note we aim at a more complete description of the thermodynamic limit 
of the Lipkin model. 



THE SPECTRUM 



A closer look at the essential features of the spectrum for large may elucidate state- 
ments made in the previous section and provides clues of what can be expected in the limit 
N oo. In Fig.l the spectra are drawn for a few values of A. As this is done for large 
values of A^ we use the 'continuous' variable x = 2k/N ranging from zero to unity and plot 
e(x) = 2Ek/N versus x. Using this scale we note the following features; for A > 1 each 
spectrum has a singularity when it crosses the critical line e = — 1. In fact the singularity 
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FIG. 1: Spectra for A = (solid line), A = 1 (dotted line), A = 5 (dashed line) and A = 10 (solid 
line). 

is a point of inflection with a zero derivative and an inflnite second derivative as shown in 
Fig. 2 and in accordance with Eq.© (for A = 1 the spectrum just touches the critical line 
with zero flrst and inflnite second derivative). In other words, the level density is singular 
at the critical line. Note that the singularity occurring at the critical line is common to all 
spectra for A > 1 irrespective of the particular value for A. For A < 1 the spectra, i.e. their 
level densities, are smooth. 

The range > e > — 1 is the normal phase while e < — 1 corresponds to the deformed 
phase. For better illustration we do not display the lower end of Fig.l, it suffices to note 

n 

that [2] for A > 1 the ground state assumes the value 



The critical line is crossed by each level sequence at a certain Xc(A) associated with A > 1; the 
larger A the larger Xc, i.e. the phase transition occurs higher up in the spectrum and higher 
above the ground state given in Eq.®. The critical line is special not only as it signals the 



<0) = --(A + -). 



(6) 
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FIG. 2: Spectra (top row) and their derivatives (bottom row) for A = 1 (left column) and A = 10 
(right column). The top row is displayed for convenience, it shows the same respective curves 
including the critical line as those in Fig.l. 



phase transition through the singularity, but it is also associated with a specific state vector 
being identical to the unperturbed ground state once the limit — > oo is attained. In fact, 
it is 

^^(A)|j, -j) = -j) + v^^lj, -J + 2) (7) 

where the second term of the right hand side vanishes in the limit. The right hand side also 
becomes independent of A when N ^ oo. Recall that for finite N all levels and state vectors 
are analytic functions of A and analytically connected in the A-plane. As a special result we 
further note that this state is perfectly localised; all other states are more or less extended 
superpositions of vectors \j,m), —j<Tn<j; this localisation was noted also in [3]. 

The singularity discussed above for the level density is found as well when a particular 
level (fixed value of x) is plotted versus A. The same pattern prevails: there is a point of 
inflection at A = Ac (associated with e{xc) = —1) with a vanishing derivative and an infinite 
second derivative in A. 
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FIG. 3: Exceptional points in the A-plane for N = 8 (diamonds), = 16 (triangles) and = 32 
(dots). The first twelve points along the real axis and the innermost arc running from the real to 
the imaginary axis are indicated by asterisks for = 96. 



EXCEPTIONAL POINTS 

The whole pattern, in particular the singular behaviour at the critical line, i.e. at the 
phase transition, can be understood from the global analytical structure of the spectrum 
as a function of the complex parameter A. The exceptional points (EP) ^| are the only 
singularities of the spectrum for the underlying model. For a hermitian problem they can 
occur only in the complex A-plane. They are square root branch points where, generically, 
two energies are connected by a branch point. For the Hamiltonian (0) all energies Ek{X) 
are analytically connected, in fact, each Ek{X) represents the values on a particular Riemann 



sheet of one analytic function 12 1 . 



In Fig. 3 we illustrate the exceptional points for a few small values of N. Two aspects are 
of importance: (i) for increasing some EPs accumulate along the real axis A > 1 and, (ii) 
at the EPs of the sequence along the real axis, starting at A = 1, the levels coalescing are 
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El and E2, E2 and £"3, E^ and E4 and so forth. 

The second point clarifies immediately some observations made in the previous section. 
Take, for instance, A = 2. From Fig.3 we read off that there are N/8 EPs along the real 
axis between A = 1 and A = 2. This implies that the Riemann sheet (energy level) with 
k = N/8, that is for x = 0.25, has a singularity at A = 2. In general terms, there is just one 
Xc associated with a particular value of A. 

With regard to the first point we note two aspects that are, however, related. In the 
limit N ^ 00 there is a dense set of EPs, say at A = 2, for x = 0.25. The singularity, being 
an accumulation of square root branch points, cannot be expected to be of a square root 
nature. The more likely behaviour is that of a logarithmic singularity. This is supported by 
a similar, albeit simpler situation where an accumulation of square root branch points 
has been demonstrated analytically to lead to a logarithmic branch point. There is further 
strong numerical evidence: a perfect fit is obtained of the curves in Figs.l and 2 using the 
trial functions 

(x — Xc)'^{ai In \x — Xc\ + a.2(ln \x — a;c|)^ + . . .) (8) 

employing three terms or less. The trial function (jSJ is used to fit the spectra of Fig.l (after 
a shift upward by one unit to make the critical line occurring at zero). In a subsequent acid 
test the derivative of (jH)) is compared with the respective curves for Ae in Fig. 2 yielding 
an excellent agreement. Note that the modulus of being used implying that the 

function is non- analytic at this point. Depending on the sign of X X Q 3j separate fit must 
be made. As it was mentioned in the previous section the same pattern is found when a 
level is considered as a function of A. A good fit is obtained replacing x and Xc by A and Ac, 
respectively, in (jH}. 

Not only do we have evidence for a very special singularity, the fact that the EPs ac- 
cumulate on the real axis indicates that the thermodynamic limit of the Hamiltonian 
cannot be a hermitian operator. Moreover, if the limit is taken either in the deformed or in 
the normal regime, the two regimes become disconnected for real values of A. 

CONCLUSION 

We summarise our findings: the thermodynamic limit ^ 00 of the Hamiltonian (^J 
leads to a specific singular behaviour associated with the phase transition. It is immaterial 
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as to whether a specific level is considered as a function of A or whether the whole spectrum 
is considered for a fixed value of A. In either case, a logarithmic singularity is found at a 
particular value of A being uniquely related to a specific level. Note that the singularity is 
non- analytic in that, on the real A-axis, the behaviour on the left hand side of the singularity 
is not the continuation of that on the right hand side. The singularity produces an infinity 
in the respective second derivative. As such, it is no longer a second order phase transition 
as it is found for finite A^. 

Furthermore, the fact that the EPs accumulate on the real A-axis in the limit, renders 
the operator nonhermitian in the limit. In fact, depending on how the limit is taken, a 
different result emerges: one operator describing the normal phase and the other describing 
the deformed phase. In contrast to the case of finite A^, the two are dissociated. This 
is in contrast to finite A^, where all energies and state vectors are analytically connected. 
In particular, when the limit A^ oo is attained, the state vector \j,—j) (see Eq.([7j)) is 
common to all spectra, irrespective of A. As a particular feature, this vector has, for A > 0, 
a non-zero overlap with any other vector of either the normal or the deformed sector. This 
is in contrast to finite A^, where yields the usual complete orthonormal basis for any 
given A. 

It ,s worth to ,ne,.,o. that the se,„,c,as*a, approach Q ha. expHct.y g,ve„ nse to 

ambiguities of a hermitian thermodynamic limit. An unambiguous limit could be obtained 
only to lowest order in ^ ~ 1/A^. 

We refrain from speculating about physical consequences of our findings. Rather we 
feel that the thermodynamic limit of the Lipkin model is unsuitable to describe a physical 
system despite its beauty for finite yet large values of A^. We see the value of this analysis, 
apart from its own sake, in a specific contribution of the general understanding of models 
describing quantum phase transitions. 

Acknowledgment 

The author gratefully acknowledges helpful discussions with Hendrik Geyer and Frederic 
Scholtz. 



8 



[1] H.J. Lipkin, N. Meshkov and N. Click, Nucl. Phys. A 62, 188 (1965); ibid. 62, 199 (1965), 
ibid. 62, 211 (1965) 

[2] P. Ring and P. Schuck, The Nuclear Many Body Problem, Springer, New York (1980) 
[3] R. Botet, R. Jullien and P. Pfeuty, Phys. Rev. Lett. 49, 478 (1982) 
[4] J.I. Cirac, M. Lewenstein, K. M0lmer, and P. Zoller, Phys. Rev. A 57, 1208 (1998) 
[5] J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, 022107 (2004) 

[6] F.G. Scholtz, B.H. Bartlett and H.B. Geyer, Phys. Rev. Lett. 91, 80602 (2004); J.N. Kriel, 

A. Morozov and F.G. Scholtz, J. Phys. A 38, 205 (2005) 
[7] F. Leyvraz and W.D. Heiss, Phys.Rev.Lett. 95, 050402 (2005) 
[8] S. Dusuel and J. Vidal, Phys. Rev. Lett. 93, 237204 (2004) 

[9] W.D. Heiss, F.G. Scholtz and H.B. Geyer J.Phys. A: Mathematical & General, 38, 1843 
(2005) 

[10] M.Z.L Gering and W.D. Heiss, Phys. Rev. C 29, 1113 (1984); W.D. Heiss and R.H. Lemmer, 

Phys. Rev. C 29, 2368 (1984) 
[11] T. Kato, Perturbation theory of linear operators (Springer, Berlin, 1966) 
[12] W.D. Heiss and A.L. Sannino, Phys. Rev. A 43, 4159 (1991) 
[13] W.D. Heiss, M. Miiller and I. Rotter, Phys.Rev. E58, 2894 (1998) 



9 



